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Abstract - Designing and optimizing cost functions and energy landscapes is a problem en- 
countered in many fields of science and engineering. These landscapes and cost functions can 
be embedded and annealed in experimentally controllable spin Hamiltonians. Using an approach 
based on group theory and symmetries, we examine the embedding of Boolean logic gates into 
the ground state subspace of such spin systems. We describe parameterized families of diagonal 
Hamiltonians and symmetry operations which preserve the ground state subspace encoding the 
truth tables of Boolean formulas. The ground state embeddings of adder circuits are used to 
illustrate how gates are combined and simplified using symmetry. Our work is relevant for exper- 
imental demonstrations of ground state embeddings found in both classical optimization as well 
as adiabatic quantum optimization. 



The embedding of energy landscapes into the ground 
state subspace of spin systems is a task commonly en- 
countered in both classical [1-3] and quantum optimiza- 
tion [4-6]. Finding a state in this subspace is equivalent 
to a wide variety of NP-complete decision problems and 
NP-hard optimization problems [1,2.7-10] which have re- 
' ceived renewed interest in the wake of adiabatic quan- 
tum computation [4, 11-15] and its experimental realiza- 
tions [16-18]. Recent works have focused on embedding 
cost functions into the ground state subspace of spin sys- 
tems [5,9,11-15,19-25] and cellular automata [21,26-28]. 
While the emphasis and techniques used in previous work 
varies, many of the fundamental results overlap. 

In this letter, we use symmetries of Boolean functions 
to unify and extend various constructions of Hamiltonians 
embedding Boolean functions into their ground state sub- 
spaces. We perform a systematic analysis of the Hamilto- 
nians embedding all two-input, one-output gates using our 
group theoretic approach. We also report on a new family 
of Hamiltonians embedding the universal logic gate NAND 
and present a new XOR Hamiltonian embedding which en- 
compass several previous results [19,20,25]. Both of our 



constructions have three free parameters providing previ- 
ously ignored degrees of freedom which could be useful 
when considering experimental constraints. Extensions of 
our symmetry arguments to larger Boolean functions are 
demonstrated using adder circuits of increasing complex- 
ity. 

While we focus on embedding circuits into the ground 
state, the application of symmetry arguments is quite gen- 
eral and can be used in the construction of Hamiltoni- 
ans for other embedding problems recently studied in adi- 
abatic quantum computing such as lattice protein fold- 
ing [24,29], adiabatic quantum simulation [30], machine 
learning [23], or search engine rankings [22]. 

Throughout this letter, we use diagonal Hamiltonians 
of N spins 
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with a = a z defined by a = |0)(0| - |1)(1|. Since the 
eigenvalues of a are ±1, we identify Boolean variable, x G 
{0, 1}, with (1 — er)/2 instead of a itself. The subscript of 
each a indicates which spin the operator acts on. Terms 
such as (jjcr j are understood as the tensor product <n ® aj. 

Limiting the Hamiltonian in eq. (1) to two-spin inter- 
actions yields the experimentally relevant [16-18] tunable 
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Ising Hamiltonian which will be our primary focus. 

The idea of ground state spin logic is to embed 
Boolean functions, / : {0,1}" -» {0, l} m , into the 
ground state subspacc, £(ify( x )), of spin Hamiltonian 
H ■■■ , <7 k ) acting on the spins a. t , a } , . . . , a k . 
As an example, consider the universal NAND gate defined 
by NAND(x,?/) = x\l y. The corresponding Hamiltonian, 
H x\/ yi& i , 03) , should have the following ground state 
subspacc 

C{H x yy) = span{|a;)|y)|sVy)} (2) 
= span{|001), |011), 1 101), |110)} 

Using the a matrices, such a Hamiltonian is given in [21] 
as, 

H xV y{cr 1 ,a2,cr 3 ) = 21 + (1 + ai + cr 2 - aia 2 ) cr 3 (3) 

This construction uses a three-spin interaction which can 
be replaced using the same number of spins and only two- 
spin interactions. This was done in [19, 20] by penalizing 
and rewarding certain interactions such that the ground 
state subspace is not altered while the higher energy eigen- 
states are. 

Now we introduce the first result of our paper: a three- 
parameter family of Hamiltonians that generalizes the for- 
mulas found in [19, 20, 25] and elsewhere. Using coeffi- 
cients labeled as in eq. (1), the constraint that one eigen- 
subspace is four-fold degenerate and contains states |001), 
1 1 1 ) , 1 101), and 1 1 10) leads to the following three equali- 
ties: 

C3 = ci + c 2 (4) 
C13 = C12 + ci (5) 
C23 = C12 + c 2 (6) 

After enforcing these constraints, the energies are 

Edegen = —C\ ~ C 2 — C 12 (7) 
£000 = 3(c 1 +C 2 + Ci 2 ) (8) 

E w = 3ci - c 2 - C12 (9) 
£100 = 3c 2 - ci - C12 (10) 
Em = 3c i2 -ci-c 2 (11) 

For ci, c 2 , and ci 2 greater than zero, the degenerate space 
is always the ground state. In closed form the three- 
parameter family of Hamiltonians encoding NAND in the 
ground state is 

Hxvy((?i, 0-2,0-3) = (ci<7i + c 2 er 2 )(l + 03) (12) 

+ (C1 + C 2 )(T 3 + C12 ^ (TiCJj 

i<j 

with ci,c 2 ,ci 2 > 0. The freedom to select these parame- 
ters could be desirable as it reduces the constraints placed 
on an experimental realization. 

The ground state energy of the NAND Hamiltonian, is 
— (ci + c 2 + C12) instead of zero. Some authors choose to 
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Fig. 1: (Color online) Ground state embedding of the half adder 
circuit, a) The half adder is implemented with a XOR gate and 
an AND gate, b) The XOR and AND gates have been substi- 
tuted by the corresponding all-NAND circuits, c) The same 
circuit has been rewritten without the redundant gates and 
labeled wires, d) Here the circuit is mapped to a network of 
seven spins, each corresponding to the seven wires of the cir- 
cuit. The thickness of each link is proportional to the two-spin 
interaction strength, while the size of each node is proportional 
to the local field strength in the two-local reduction. The pa- 
rameters used for the NAND gate Hamiltonian given in eq. (12) 
are ci = c- 2 — C12 = 1. 



consider positive semi-definite Hamiltonians, however the 
addition of multiples of the identity does not alter energy 
differences within the landscape of the problem and we 
choose not to enforce this constraint. 

As the NAND gate is universal for the construction of 
logic circuits, the NP-complete problem CIRCUIT-SAT, 
where the question "Is there an input corresponding to the 
output of logical one?" is embedded using only positive 
couplings and positive local fields. This leads to an alter- 
native proof that finding the ground state of spin Hamil- 
tonians with anti-ferromagnetic couplings in a magnetic 
field is NP-hard [2]. 

Let us turn to an illustration that shows how to use 
the Hamiltonian in eq. (12) to construct more complex 
functions. Naively, it may seem a separate spin must be 
included for each wire originating from a FANOUT oper- 
ation [19-21]. However, this is not the case; instead the 
same spin may be used for the input to as many gates as 
desired. As an example, in fig. 1, an all-NAND half adder 
circuit is converted to a spin Hamiltonian using eq. (12). 
We will return to this example at the end of the letter as 
an application of our symmetry considerations. 

An important consideration for this model is the in- 
put and output of the circuit. To extract data from this 
system, single spin projective measurements can be used. 
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Inputs are set using an additional Hamiltonian 



F U F 2 ,R 



inputs 



Hi, 



(13) 



which forces the fc-th bit to take the value x k S {0, 1}. 

There are certain symmetries of Boolean functions from 
which we can infer properties of the class of Hamiltonians 
that have the Boolean function embedded in the ground 
state subspace. 

To limit the scope of our initial discussions, we will re- 
strict our attention to Hamiltonians containing only two- 
spin interactions and to the set of the 16 two-input, one- 
output gates. 

Each of the two-input, one-output gates is dchned by 
its truth tabic: 
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with bi £ {0, 1}. There are 16 choices for the vector b = 
[bi,b2,b3,bi]. The corresponding Hamiltonian, Hb, must 
have ground state subspace 

C(H b ) = {|006i>, |016a>, |106s>, |11&4» (14) 

Thus, there are 16 relevant ground state subspaces, each 
corresponding to one of the truth tables. 

The symmetry operations on truth tables must treat the 
output bit differently in order to remain in the space of 
the 16 truth tables. Thus, we consider (i) bit flips of any 
of the spins and (ii) swaps of the two inputs giving the 
following symmetries: {e, Fi, F 2 , F 3 , R\ 2 }- Here e is the 
identity operation, Fi is the spin- flip operation (negate), 
and R12 is the spin-swap operator (permute). The action 
of the latter two operations on spins is defined via 



F t oa 3 = (1 



25ij)(jj 



(15) 

Rij ° Cfc = o-j5 k i + o-,4j + Cfc(l - hi - hj) (16) 

The group G can be presented as 

G = (Ri2,Fi,F 3 ) 



(17) 

where (•) indicates a set of generators. Defining relations 
of the group arc R\ 2 = F? = F% = e, R 12 Fi = F 2 Ri 2 , 
FiFs = F 3 Fi and R12F3 = F 3 Ri 2 . From these relations, 
or alternatively from the cycle graph, the group is of order 
16 and is isomorphic to D4 x Z 2 , where D4 is the symmetry 
group of the square and Z 2 is the cyclic group of order 2. 

The action of G on the set of 16 truth tables is depicted 
in fig. 2. Four orbits are found under action of the group: 

{0,1}, 
{x,y,x,y}, 

{x V y, x V y, x V y, ■ ■ ■ , x A y} 
{x®y,x== y}. 
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Fig. 2: (Color online) The action of D4 x Z2 on the 16 Hamil- 
tonians corresponding to truth tables of two-input, one out- 
put functions. The Hamiltonians can be converted to any 
other Hamiltonian in the same orbit by applying the spin-flip 
(negate) Fi or input-swap (permute) R12 operations. The sym- 
metry operations that leave the ground state subspaces of each 
Hamiltonian invariant (the stabilizer subgroup) is written on 
the perimeter of each rectangular region. Orbits a), b), c), d) 
are explained separately in the text. Each of these orbits re- 
quires an additional spin for a Hamiltonian embedding using 
only two-spin interactions: orbit a) requires a single spin, b) 
two spins, c) three spins, and orbit d) requires four spins. 



These classes are depicted in fig. 2a, 2b, 2c, 2d, respec- 
tively. These classes correspond to different NPN (negate- 
pcrmutc-ncgatc) classes [31,32]. Interestingly, each orbit 
requires a different number of spins to implement when 
considering only two-spin interactions. We examine each 
in turn. 

First, consider the constant functions with bi — c and 
c G {0, 1}. Since these functions do not depend on x nor y, 
there is no need to couple either to the third spin. Hence, 
the Hamiltonian in eq. (13) can be used. According to the 
group action depicted in fig. 2a, given the Hamiltonian for 
Hzero corresponding to bi = 0, the action of F3 transforms 

Hze.ro to H one . 

Second, for each of the functions, bi — x, bi — y, bi = x 
and bi = y 7 the output bit only depends on one of the two 
inputs. The other input is extraneous, so the gate only 
requires two spins to implement. The truth tables can be 
embedded using variations of the COPY gate previously 
introduced in [19-21]. The general fc-COPY gate forces k 
bits to take the same value and the corresponding diagonal 
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operator 

Hk-copy = - - a i a i (18) 

acting on fc-spins possesses a ground state subspacc 

£(fffc-COPY) = span{|0)® fc , |l)® fe }. (19) 

If we are concerned with constructing a Hamiltonian us- 
ing a physical set of spins, the spatial locality could play 
an important role as coupling of distant spins may not be 
possible. In this case, the fc-COPY gate could be useful for 
spatially distributing intermediate results of the computa- 
tion. The action of F 4 or F3 transforms H x = H 2 -copy( x ,z) 
into the Hamiltonian H x , as shown in fig. 2b. 

The third class of functions to be considered is x V y, 
x A y and all possible negations of the two inputs. Our 
general formula for x V y is given in eq. (12) and using the 
symmetry operations from group G, see fig. 2c, all other 
gates in this orbit can be derived using three spins with 
two-spin interaction terms (see the appendix for additional 
formulations) . 

The last orbit of functions, XOR (x © y) and its logical 
negation EQUIV (x == y), cannot be embedded in the 
ground state subspace of a three spin system using only 
two-spin interactions; it requires a fourth ancilla spin to 
implement using only pairwise interactions. If restricted 
to three spins, the gate XOR (©) requires a three-spin 
interaction. 

H xS ,y(ai,(72,cr 3 ) = -<7ia 2 <T 3 (20) 

The inability to create this operator acting on three 
spins with two-spin interactions can be demonstrated alge- 
braically or graphically using Karnaugh maps [19,25]. For 
XOR, the stabilizer subgroup is generated by FiFj and 
i?i2, see fig. 2d. When considering the ancilla spin, 04, 
there is an additional F4 symmetry that leaves the truth 
table unchanged. 

Beginning with the swap-symmetric operators M z = 
J2i a i an d M zz = Yli<j c r i° T 3 '■> we wr ite the most general 
swap-symmetric Hamiltonian over four spins restricted to 
two-spin interactions as 

H R = r z M z + r zz M zz + a 4 (r 4 + r z4 M z ). (21) 

Suppose that the coefficient vector R = [r z , r zz , r 4 , r z4 ] 
gives a valid XOR Hamiltonian. Then we can act with F 4 
to get a second Hamiltonian that also preserves the ground 
state subspace with coefficients R' = [r z , r zz , — r 4 , — r z4 \. 
In references [19] and [20], this F 4 symmetry connects 
the decompositions given as R = [1,-1,-2,2] and R = 
[1,-1,2, —2] in the respective papers. Furthermore, since 
the ground state subspace is symmetric with respect to 
FiFj, there are an additional six Hamiltonians with log- 
ically equivalent ground state subspaces. For example, 
beginning with H x ^ y corresponding to R = [1, —1, —2,2] 



z = f(x,y) H f{x ^ y) (a 1 ,a 2 ,cr 3 ,a 4 ) 

Constant functions 

2=0 H zero = (1 - o- 3 ) 

Copy-type functions 

Z = X H x = (1 - (71 cr 3 ) 

AND, OR, ... , NAND, NOR functions 

_ w _ H xV y= (ciCTi + c 2 er 2 )(l + cr 3 ) 
z = xVy 1 / 1 \ 1 ^3 

+ (ci + C 2 )cr 3 + C12 }^ l<j ViVj 

XOR and EQUIV functions 

ffi H x@y = H xA y(ai,a 2 ,cr 4 ) - a 3 

z — x ty y r, 

+ 0103 + CT20-3 + 2(73(74 



Table 1: Summary of representative Hamiltonians from each 
orbit under the action of the symmetry group. Spin one and 
two correspond to the two inputs while spin three corresponds 
to the output. The fourth spin is an ancilla spin needed only 
for the implementation of XOR and EQUIV. In the AND, OR, 
. . . , NAND, NOR family, the sign of the coefficients determines 
which gate on this NPN orbit one obtains, as detailed in the 
appendix. We have only shown four Hamiltonians and the 
remaining 12 Hamiltonians as well as additional Hamiltonians 
with different excited states are related via the action of the 
group D 4 x Z2 as depicted in fig. 2. 

and using symmetry operation F\F 2 results in 

F 1 F 2 oH x(Sy = 2cr 4 (-(7i -CT2 + CT3) (22) 

+(7l + (72 — CT3 — 2(74 
+ (<71<72 — (72<7 3 - (7i(73) 

with the same ground state subspace. Note that this 
Hamiltonian is not of the same form of eq. (21) like those 
given in [19,20]. 

To extend the XOR Hamiltonians previously listed to 
a parameterized family of Hamiltonians, we rearrange 
eq. (21) with R = [1, -1, -2, 2] as 

H x(Sy = -((71 +<7 2 )(1-<T 4 ) (23) 

— 2(7 4 + ((7l<7 2 + 0-1(74 + O2O-4) 

— CT3 + CT1O-3 + (72(73 + 20-30-4 

Comparing with eq. (12) and using fig. 2c, we can sim- 
plify this equation using H xA y(a 4 ,a 2 ,a 4 ) = F 4 F 2 F 4 o 
H x \zy{ai,a 2 ,a 4 ) evaluated at ci = C2 = C12 = 1- Gen- 
eralizing to other values of ci,C2, and C12, we arrive at 
the following three-parameter family that preserves the 
ground state subspace of XOR 

F[ x ® y = H. xA% y((j 1 , 0-2, 04) - 0-3 (24) 
+O10-3 + 020-3 + 20-30-4 

By examining the excited state structure of eq. (23), we 
find that in the parameterization of H x/ \y the coefficients, 
ci, c 2 , C12, must be greater than 1/2 instead of strictly pos- 
itive. 

Our work has direct relevance to recent experimental 
realizations of adiabatic quantum computing in supercon- 
ducting qubits [17,18] and ion traps [16] where controllable 
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Fig. 3: (Color online) The half adder spin Hamiltonian that 
arises from the four spin decomposition of the XOR Hamilto- 
nian which simplifies the construction from fig. Id. Dashed 
links represent negative interactions and checkerboard shad- 
ing indicates a negative local field. The size of the nodes 
and the thickness of the edges are proportional to the fields 
and interaction strength (on spins A and B there is no local 
field). The parameters used for the AND gate and XOR gate 
are c\= C2= c\% = f . 

couplings between spins can be used to embed problems 
into the target Hamiltonian of the evolution. Since both of 
these experimental systems are limited to two-spin inter- 
actions, our decomposition for XOR provides an effective 
three-spin interaction which is experimentally realizable. 

In table 1, we summarize our results for Hamiltonian 
embeddings of two-input, one-output Boolean functions. 
While we have restricted attention to diagonal Hamilto- 
nians, future work could consider transformations where 
the ground state is preserved but the Hamiltonian obtains 
off-diagonal elements. 

Now we return to the half adder example from fig. 1. 
With our constructions, we can directly implement it using 
the XOR and AND gates, 

Hha = H xS)y (a a, &b , cr a , as) + H xhy (aA, o Bl ac)- (25) 

Here a a and <tb correspond to the inputs to be summed, 
a a corresponds to the XOR ancilla bit, and as and ac 
correspond to the sum and carry bits. As depicted in fig. 3, 
the new spin Hamiltonian uses two less ancilla spins than 
our earlier construction and now has six free parameters. 
Additional degrees of freedom arise from the D± stabilizer 
subgroup of the XOR Hamiltonian and the Zi stabilizer 
subgroup of the AND Hamiltonian. 

The symmetry group of Hha can be inferred from the 
symmetries of the component Hamiltonians using a di- 
rect product structure. For a general circuit Hamilto- 
nian composed of gate Hamiltonians acting on subsets of 
spins, H = Hi, the stabilizer subgroup is the direct 
product of the stabilizers for each of the Hamiltonians in 
the sum. The direct product group action is defined as 
(<7i, (72, • • ■ 9n) ° H — ^2 gi o Hi. If g is in the intersection 
of all stabilizer groups (the diagonal subgroup), then goH 
will have the same ground state subspacc as H . 

Additional symmetries arise after partitioning the bits 
into output and ancilla bits. We can expect the symme- 
tries of the Boolean function being embedded to be pos- 



sessed by the resulting Hamiltonian. However, the symme- 
try group composed of the gate-local symmetries preserves 
the full ground state subspace including the values of the 
ancilla bits. The symmetries of the Boolean function be- 
fore being decomposed into logic gates will arise as global 
symmetries that cannot be obtained from the gate-local 
symmetries of the individual gates. For instance, if a a 
corresponds to an ancilla spin, then inverting this bit in 
each circuit component leaves the ground state subspace 
invariant. That is, H and (F a , F a , ■ ■ ■ , F a ) o H embed the 
same Boolean function. 

As a further illustration of the distinction between 
global and gate-local symmetries, consider the full adder 
corresponding to a Boolean function which adds binary 
summands A, B, and carry- in bit Ci n . The permutation 
of the input bits and the carry-in bit is a symmetry of the 
full adder Boolean function. However, such a permuta- 
tion is not a gate-local symmetry of the sub-Hamiltonians 
used in the circuit embedding, see the appendix for de- 
tails. This is because the values of the ancilla spin within 
the ground state subspace is not preserved under this per- 
mutation. Thus, the local symmetries do not determine 
all possible symmetries when some bits are considered as 
ancillas. 

As a final example of ground state spin logic, fig. 4 shows 
the spin Hamiltonian of the ripple carry adder for four-bit 
binary numbers. The figure shows the network for both 
an implementation with only NAND gates in fig. 4a and an 
implementation with XOR, AND, and OR gates in fig. 4b. 
The second construction allows a decrease in the number 
of ancilla spins and provides 51 free parameters. Addi- 
tionally, as shown in the appendix, the symmetry group 
of the second implementation has at least 2 31 elements. 
Another salient feature is that the average degree of the 
spins changes from 3.85 in the all-NAND case to 4.22 in the 
second implementation. Explicitly listing the free param- 
eters and the symmetries that preserve the ground state 
subspace is an illustration of how our approach gives ex- 
perimentalists and theorists systematic methods to find 
additional degrees of freedom. 

An important step towards large scale experimental re- 
alizations of the techniques presented in this paper will be 
the adiabatic implementation and characterization of the 
elementary logic gates. In the case of XOR, this Hamilto- 
nian will allow one to realize an effective three-spin inter- 
action by using only two-spin interactions and introducing 
an ancilla spin. Such an interesting example is in line with 
current experimental capabilities [16-18]. 

* * * 

The authors would like to thank V. Bcrgholm and Z. 
Zimboras for helpful discussions and M. Allcgra and J. 
Roland for carefully reading the manuscript. JDW ac- 
knowledges support from NSF (No. 1017244) and thanks 
the Visitor Program at the Max-Planck Institute for the 
Physics of Complex Systems, Dresden where parts of this 
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'cout 



b) 



Fig. 4: (Color online) Ripple carry adder. The figure shows the network of spins corresponding to a ripple carry adder with 
four bits. The ripple carry adder is composed by one half adder and three full adders; in yellow it shows the input spins from 
the four bits binary numbers A = and B — Ylt—i Bi2 % ; while the sum spins, Si are drawn in purple. Carry bits are 

labeled as d. The direction of the sum is from left to right. Fig. a) shows a ripple carry constructed with only NAND gates 
and parameters ci = C2 = c\% = 1, while b) shows the same adders built with XOR, AND and OR gates. 



work were completed. 
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Appendix. 

Hamiltonians embedding full adders. Wc provide the 
characterization of the full adder [33] necessary to con- 
struct the ripple carry adder shown in the main text. This 
affords us an opportunity to explore the network proper- 
ties of the adders circuit family with well known construc- 
tions and optimized solutions [34,35]. 





a) Full adder circuit 




b) Full adder circuit with NANDs 



dcg. centr. + 
shortest path centr. X 



+ + + 

X 



+ + 



+ + 
X + + 



XX 



s, X + + 

X X 



c) Full adder network 



X-X 

ABC m 5 6 7 SC 
bit 

d) Centrality of nodes 



Fig. 5: (Color online) a) Full adder circuit, b) Transcription 
with only NAND gates, c) The spin network representing the 
all-NAND circuit of the full adder. The input spins A, B and 
the lower carry-in spin C are depicted in yellow while the out- 
put spins of sum (S) and carry-out (Cout) are colored in purple. 
The green nodes represent ancilla spins, d) Degree centrality 
and shortest path centrality of network nodes. Notice that the 
input and output spins (denoted A = 0, B = 1, d n = 2 and 
S = 12, C ou t = 13 respectively) are the least central while the 
most central nodes are the ancilla spins five, six and seven. 

In order to sum arbitrarily large binary numbers, the 
half adder circuit needs to implement the bit carrying op- 
eration. The full adder circuit introduces this operation 
with a third input bit, accounting for the lower level carry 
bit Cm- Fig. 5 shows the network associated to this cir- 
cuit, where the inputs bits A, B and Cj„ are in yellow and 
the output bits S and C ou t are in purple. The network 
corresponds to a circuit with only NAND gates. The two- 
spin interactions and the local fields are then all positive 
valued. From eq. (12), we have three free parameters for 
each NAND gate, giving 3 • 5 = 15 free parameters. Each 
NAND Hamiltonian is also symmetric under the action of 



the symmetry group {e, -R12}, giving a symmetry group for 
the whole Hamiltonian of at least 2 5 elements. The Hamil- 
tonian uses nine ancilla spins to build the truth table of 
the full adder, resulting in nine new symmetries, labelled 
in the main text as {F a : a labels an ancilla spin}. The 
action of the latter changes the ground state subspace but 
the resulting system still describes the original problem. 




out 



Fig. 6: (Color online) Full adder spin network for the circuit 
in fig. 5a. This construction reduces the number of ancilla 
spins from nine in the case in fig. 5d to five. In this case, the 
interactions are not all of the same sign. 

The number of ancilla spins can decrease using the stan- 
dard XOR, AND and OR gates of fig. 5a. Fig. 6 shows the 
spin network associated to this circuit. This Hamiltonian 
presents at least 2 9 symmetries arising from the single-gate 
symmetries. 

To enhance comprehension of the spin Hamiltonians, 
we compute some well known complex networks mea- 
sures [34] . The node centrality of the resulting network, in 
fig. 5d, suggests that the input and outputs spins are the 
least central both for local (degree centrality) and global 
centrality measures (shortest path centrality). Here the 
degree centrality D k of node k is denned as: 



N - 1' 



(26) 



where dk is the degree of node k and TV is the number 
of nodes of the network. The shortest path centrality is 
denned as: 

SP ikj 



1,3 



SPi. 



(27) 



where SP%j represents the number of shortest path be- 
tween nodes i and j and SPikj is the number of those 
paths passing through node k. Nodes with higher cen- 
trality can be thought as the network bottleneck between 
input and output spins. 

Additional calculations for the ripple carry adder. The 
sum of two binary n-bit numbers x and y, can be carried 
out by concatenating n full adder circuits yielding the rip- 
ple carry adder. This circuit implements a cascade: the 
carry-out bit of each full adder will be used as the carry-in 
bit for the next one. The first full adder has logical zero 



p-7 



J. D. Whitfield et al. 



as the carry-in bit. Alternatively, it can be completely re- 
placed by the half adder circuit. Fig. 4 of the main text, 
shows the network associated to a four bit ripple carry 
adder used to sum two four-bit binary numbers. In fig. 4a 
the Hamiltonian is built from NAND gates, providing all 
non-negative local fields and interactions and resulting in 
46 spins and 86 links. The starting circuit contains 38 
NANDs, each of them are associated to a Hamiltonian, see 
eq. (12), which depends on three free parameters. Thus, 
the Hamiltonian of the whole circuit has 3 • 38 = 114 free 
parameters. As in fig. 2, the stabilizer subgroup of the 
NAND gate contains only two elements: 

stab(NAND) = (R 12 ) ~ Z 2 (28) 

and generates a group of symmetries for the full Hamil- 
tonian of at least 2 38 elements. Fig. 4b shows the same 
circuit built using seven XOR, seven AND and three OR 
gates. This implementation yields a network with only 32 
bits and 65 links. The gate Hamiltonians have three free 
parameters each for a total of 3 ■ 17 = 51 free parameters. 
In this case, the stabilizer subgroup of both OR and AND 
is also generated by i?i2, while for the XOR gate we have: 

stab(XOR) = (F 1 F 2 ,F 1 F 3 ,R 12 ) ~ D 4 (29) 

with eight elements. Thus, total symmetry group of the 
Hamiltonian contains at least 2 10 ■ 8 7 = 2 31 elements. 

Fig. 7 shows the centrality of each spin of the net- 
works in fig. 4. We note that, in the second implemen- 
tation of the ripple carry adder, the resulting network is 
slightly more connected. The average degree centrality 
((Di) — 0.131) in this case is higher than in the implemen- 
tation with only NANDs ((A) = 0.083). The variance is 
also higher in the second example, var(D^) = 0.041, than 
in the first, var(A) = 0.021. In both implementations, the 
most important spins as identified by the global measure 
of shortest path centrality are the spins on the backbone 
of the circuit. In particular, the carry bits, Ci, have high 
centrality as they connect subnetworks which would oth- 
erwise be disconnected. 

General formulas for orbit of NAND under -D4 x Z 2 . 
From the main text, the Hamiltonian for NAND is 

H = (ciCTi + c 2 er 2 )(l + 03) + (ci + C 2 )<7 3 + c 12 M zz 

This can be written as 

H = C\ 2 {0\0 2 + (71(73 + 0-20-3) (30) 

+0-3(01(1 + 0-2(1 + 0-3))) 
+o- 3 (c 2 (l + 0-2(1 + 0-3))) 

The energy shift to ensure that the ground state is also 
the null space is 

Ci + c 2 + C12. 

We see that the symmetry in variables o~\ and a 2 breaks for 
c i 7^ c 2j yet the ground state subspace remains invariant. 



As mentioned in the main text, this degree of freedom 
could be desirable as it reduces the constraints placed on 
an experimental realization. 

Wc can identify NAND as a point on the orbit of 
the NPN class by considering three indicator variables, 
x,y,z e {0, 1}. We write 

H = ((-l)' T Cl ai+c 2 (-lF(72)(l + (-l)V3) 

+ (-1)* (Ci + C 2 )(73 + Ci 2 ((-ir + Vi(72 

+ (-lf+*<7i<7 3 + (-l) y+z or 2 <T 3 ) 

The Hamiltonian for NAND is recovered by setting 
x = y = z = 0. The rest of the orbit is picked out by 
assigning other values to x, y, z. The ground state energy 
as a function of x,y, z is given as 

E gs = Cl2 ((-ir+y - (-ir+ z - (-i)y +z ) 
-(ci + c 2 )(-iy + 
+(1 - (-rr)((-rf Cl + (_i)»oa) 

A meaningful experimental demonstration showing the 
capabilities to realize every gate in the orbit could be 
performed by realizing each of the eight Hamiltonians 
x,y, z E {0, 1}. Apart from characterizing the degenerate 
ground space of each of the eight gates, the experiments 
can also be modified slightly to correspond to instances 
of adiabatic search algorithm as follows. The output of 
each gate can be set to either logical-zero or logical-one, 
see eq. (13). Successful adiabatic annealing would then 
return the associated inputs to the circuit. 
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a) Ripple carry network with only NAND gates 
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Fig. 7: (Color online) Centrality measures for the ripple carry adders of fig. 4. Graphs a) and b) correspond to networks in 
fig. 4a and 4b respectively. In the first implementation, the number of spins used is 46 with an average degree centrality of 
0.083, while in the latter the same problem is embedded with a lower number of spins, 32, but with a higher average degree 
centrality of 0.I3I . 
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